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A.  Introduction 


Consider  a layer  of  fluid  of  conductivity  cr  , 
bounded  at  the  surfaces  z » 0 and  z = t by  semi-infinite 
insulating  media.  Let  two  parallel  bare  metal 
cylindrical  electrodes  of  radius  a,  where  a 
the  surface  z =0,  at  x = + R,  extending  from  y = 0 to 
y = L.  A low-frequency  generator  feeds  the  electrodes 
at  the  ends  y « 0 . See  Fig . 1 . 


Fig-  1 
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Now,  introduce  a relatively  small  sphere  of  radius 
rQ  and  conductivity  eg  in  the  region  between  the 
electrodes,  in  the  plane  2=0.  The  electrodes  will 
previously  have  produced  a current  field,  and  therefore 
an  electric  field,  in  this  region,  and  the  presence  of  the 
sphere  will  produce  a distortion  in  this  electric  field. 

It  Is  the  problem  of  this  report  to  calculate,  under 
suitable  restrictions,  the  resulting  change  of  electric 
potential  at  the  points  A,  B,  C,  D,  where  A,  B,  C,  D 
form  a square  array  between  the  electrodes,  in  the  plane 
z = 0. 

To  arrive  at  the  fields  by  the  complete  solution  of 
Maxwell1 s equations  for  this  special  case  would  be 
extremely  laborious.  Instead  we  shall  assume  (1)  that  the 
frequency  of  the  generator  Is  so  low  that  a direct- 
current  approximation  is  valid,  (2)  that  the  potential 
along  the  feeder  electrodes  varies  very  slowly  with  y 
so  that  the  field  in  any  transverse  3lice  y,  y + dy, 
may  be  obtained  under  the  assumption  that  the  electrodes 
are  equipotentials , (3)  that  the  electrodes  are 
infinitely  long  so  that  end  effects  are  negligible. 

With  the  above  assumptions  the  primary  field  may 
be  obtained  very  simply  from  previous  work  of 
D.  D.  Foster  (TR  12}^  ^ who  has  worked  out  the  direct- 
current  flow  across  a long  conducting  strip  due  to 


equipotential  electrodes  along  opposite  sides.  The 
field  with  the  sphere  in  place  is  then  given,  in  the 
d.c.  case,  for  any  point  outside  the  sphere,  by  the 

* 

sum  of  the  primary  field  and  the  field  due  to  a dipole 
of  a certain  strength  with  its  center  at  the  center  of 
the  sphere  (the  sphere  Imagined  removed)  and  its 
direction  parallel  to  the  direction  of  the  original 
field.  The  strength  of  the  dipole  is  simply  related 
to  the  strength  of  the  primary,  or  Inducing,  field  in 
its  neighborhood. 

Proceeding  in  this  way,  we  hope  to  obtain  an 
order  of  magnitude  calculation  for  low  frequencies. 


B.  The  Calculation 

Let  t = 10  meters 
R = 100  meters 
L =»  300  meters 
<r  » 3 mho/meter 


* Note  that,  since  the  medium  is  conducting,  there  can 
be  no  accumulation  of  charge  in  the  d.c.  case.  Thus 
we  have,  not  an  electrostatic  problem,  but  a steady- 
current  problem.  TEeHJIpoIe- may  be  regarded  either 
as  an  electric  dipole  or  as  a current  dipole,  the 
two  being  related  through  the  conductivity  of  the 
medium.  We  suggest  that  thinking  in  terms  of  a 
current  dipole  tends  to  less  confusion. 
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Let  the  locations  of  the  points  A,  B,  C,  D be 
A:  x - +50  ns  y = 100  m 

B:  x **  -50  w y = 100  id 

C:  x = +50  ra  y - 200  m 

D:  x =*  -50  m y = 200  id 

In  TR  12  it  was  found  convenient  to  assume 
electrodes  of  semi-circular  cross-section.  The 
electrodes  in  the  actual  case  to  which  vie  wish  to 
approximate  are  0000  gauge  copper  cables  with  a roughly 
circular  cross-section  of  outside  diameter  1.33  cm.  We 
shall  assume  for  simplicity  that  the  best  choice  of 
radius  for  the  semi-cylindrical  electrode  to  give  a 
reasonably  equivalent  spreading  resistance  to  the 
actual  case  is  that  radius  which  will  give  an  equal  ex- 
posed periphery  in  the  conducting  medium.  Therefore 
we  choose  the  radius  of  the  seml-cyllndrical  electrode 
to  be 

-2 

a ® 1.33  x 10  meters 

Since  the  electrodes  are  fed  by  the  generator  at 
one  end  they  may  be  regarded  as  a parallel-wire  trans- 
mission line.  The  d.c.  transmission  line  equations  for 
the  voltage  V across  the  lines  and  the  current  I in  the 


lines  are 


*.v 


T 

(Mi 


= -£v 

3y 

Prom  handbooks  we  obtain  a value  for  the  average 
resistance  per  unit  length  of  the  transmission  line  (or 
of  two  meters  of  0000  gauge  copper  cable)  of 


$ = 3. 2 x 10“4 


ohm/raeter 


From  TR  12,  p.  10,  we  obtain  the  shunt  conductance  per 

« 

unit  length  of  the  transmission  line 


0 728  tnho/mei'er' 


If  we  neglect  any  correction  for  the  termination 
of  the  line  at  y = L,  we  obtain  for  the  solution  of  the 
transmission  line  equations 


(4) 


Substituting  the  values  for  and  <0/  , 

-0,0064  y 

V = e 

~~V~ 

o 


Values  of  V/Vq  for  various  values  of  y are  shown  in  column  2, 
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Fig.  2.  It  will  be  seen  that  the  value  of  V/V  at 

o 

y a L is  not  entirely  negligible,  being  about  1/7  of 
the  input  voltage,  and  therefore  that  an  appreciable 
end  effect  might  be  expected.  However,  on  account  of 
the  tendency  of  the  current  flow  lines  between  the  out 
ward  ends  of  the  electrodes  to  spread  out  into  the 
medium,  a better  approximation  to  the  terminal  con- 
ditions probably  would  be  to  assume  some  finite 
terminating  conductance  rather  than  to  assume  an  open 
end,  and  thus  the  ignoring  of  end  effects  altogether 
may  be  as  justifiable  an  approximation  to  make  as  any. 


TR  12,  p.  9,  gives,  with  a simple  substitution. 


for  the  electric  field  anywhere  in  the  plane  z = 0. 

This  formula  was  arrived  at  under  the  assumption  that  the 
electrodes  were  equipotentials . But,  proceeding  as  out- 
lined in  part  A,  we  merely  substitute  V from  equation  l1!-) 

for  V in  equation  (5). 
o 

A further  simplification  can  be  made  if  we  note  that, 
except  for  points  near  the  lines  (x  = + R),  the  square 
bracket  in  (5)  is  of  order  unity.  Thus,  if  we  restrict 
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ourselves  to  |xj  ^ 75  meters,  we  have,  wit 

negligible  error. 


E =■ 


<r  t 


(6) 


that  is,  E is  independent  of  coordinate  x,  for  jxj 

75  meters  . Substituting  numbers, 

_3  -6.4  x 10  3y 

E ==-*?-  .27  x 10  ; e 


(6a) 


The  values  of  E are  calculated  in  column  5 of  Fig.  2. 

T" 

o 

Now  it  may  be  shown  that  the  moment  of  the  current 
dipole  due  to  a perfectly  conducting  (p.c.)  sphere  of 
radius  rQ  in  a uniform  conducting  medium  is  given  by 


M. 


• c. 


4rrr  o-  E 


(7) 


and  the  moment  of  the  current  dipole  for  a perfectly 
insulating  (p.i.)  sphere  is 

Mp.u  = ~zrrrf<rE 


(8) 


The  use  of  equations  (7)  and  (8)  is  complicated  by 
the  proximity  of  the  sphere  to  the  insulating  boundary, 
the  plane  z = 0.  The  correct  procedure  may  be  made  clear 
by  the  following  argument.  Consider  a conducting  sphere 
in  a uniform  field  E,  the  field  being  in  the  x-dix-ection. 
The  resulting  field  configuration  will  be  symmetrical 
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with  respect  to  a plane  through  the  diameter*  parallel 
to  the  x-direction.  The  same  configuration  as  far  as 
a point  in  the  upper  half  space  is  concerned  will  be 
given  by  the  theory  of  images  if  one  considers  a half- 
sphere in  a field  E resting  on  a perfectly  reflecting 
plane.  That  is  to  say,  where  we  have  a reflecting 
plane,  the  correct  result  for  the  dipole  strength  in 
the  region  above  the  reflecting  plane  will  be  given  by 
formulas  (7)  or  (8)  if  we  postulate  that  the  object 
resting  on  the  plane  z = 0 is  a half-sphere  of  radius 
rQ.  When  the  object  introduced  in  the  actual  case  is 
some  other  shape  it  seems  reasonable  that  the  best 
approximation  should  be  reached  if  the  half-sphere  is 
chosen  to  have  an  equal  volume.  Thus  if  we  introduce 
in  the  physical  case  a sphere  of  diameter  one  meter 
(radius  0.5m)  we  choose  ro  in  equations  { 7)  and  (8) 
to  be 

£ = * O S 


Therefore  we  have,  for  the  physical  sphere  one  meter  in 
diameter 

M = 4TT X 2 (o-  S? O- E 

C.  ' y 


(7?0 


tt  <r  E 
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1. 

2. 

3. 

4. 

y 

-0.0064y 
y/V  **e  J 

0 

E/V  *=.  00427V 
° *0 

Mp.c.  {1  ' )*>6.85*10*4E 

vo 

0 

3 .0 

4.27X10"3 

2.67x10 

25 

0.852 

3. 64x10 

2.27 

50 

0.726 

3.10 

1.94 

75 

0.619 

2.64 

1.65 

100 

0.527 

2.25 

1.41 

125 

0.449 

1.92 

1.20 

150 

0.383 

1.63 

1.02 

175 

0.326 

1.59 

0.87 

200 

O.278 

1.19 

0.74 

225 

0.237 

1.01 

0.63 

250 

0.202 

0.86 

0.54 

275 

0.172 

0.75 

0=456 

300 

0.147 

0.63 

0.394 

Fig.  2 
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We  need  consider  only  the  perfectly  conducting  sphere, 
since  the  values  for  the  Insulating  sphere  are  simply 
half  in  magnitude  and  opposite  in  sign. 

The  potential  of  a current  dipole  of  moment  M, 

( p) 

located  at  (x*,  y* , 0)  is,  from  TR  18'  , p.  4. 


where 


(10) 


p/Q  , the  function  y/  (x,  y,  z.)  for  x3  = y3  =0  has  been 
computed  by  Mrs.  Lois  Edelstein  by  methods  developed  in 
TR  18.  For  other  values  of  z»  and  y»,  "p"  raay  be 

obtained  by  simple  translation.  Fig.  4 records  the 
values  of  ~p  at  the  points  A,  E,  G.  D,  in  the 
following  way.  A plan  view  of  the  configuration  is 
drawn,  divided  into  25-meter  squares.  The  center  of 
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each  square  marks  the  point  at  which  a dipole  has 
been  posited  for  the  computation  of  its  potential  field 
at  the  points  A,  B,  C,  D.  The  relative  values  of  the 


po  ven  v , uhclw  xq  i the  V 


Mirte  at- 


*,/ 

I /✓ 


T.r\&  i 013.1* 


points,  are  recorded  in  the  four  quadrants  of  each 
square,  as  shown  in  Pig-  3- 


marks  point 

(x'j'.O),  the  point 
of  location  of  a 
sphere  - 


Pig.  3 


It  remains  to  multiply  each  number  of  Fig.  4 by  the 

appropriate  value  of  (~l  ) to  obtain  the  values 

4TT<r^2tJ 


of  Cp  , the  potential  due  to  a conducting  or 
insulating  sphere  in  the  ccnaueting  medium.  Substituting 
from  (7 a),  for  the  perfectly  conducting  1 meter  sphere-*, 

/%.c.  (jjf  = rr<rE  f M2 
477-0-  2r/  4-rr<r  '2?' 


- 6-2S * /0~4E 


(since  t=10  m) 
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The  values  of  this,  in  terras  of  V , for  various  y, 
are  listed  in  Fig,  2 column  4.  The  values  of  0 „ 
in  terms  of  V , are  displayed  in  Fig.  5»  similarly  to 
Fig.  4. 

C.  Limits  of  the  Results 

The  potentials  at  A,  B,  C,  D due  to  the  introduction 
of  a sphere  can  vary,  as  noted,  from  the  values  in 
Fig.  5 for  a perfectly  conducting  sphere,  to  values  of 
half  the  magnitude  and  opposite  in  sign  for  a perfectly 
insulating  sphere,  with,  of  course,  the  possibility  of 
zero  values  for  a sphere  of  average  conductivity  equal 
to  that  of  the  surrounding  medium. 

Otherwise,  the  potentials  depend  only  weakly  on  the 

value  of  c r chosen  for  the  conducting  medium,  for  a 

constant  V , since  it  can  be  seen  that  G~  cancels  out 
o 

from  the  equations  except  where  it  appears  to  the  1/2 

power  in  the  exponent  in  equation  (4).  V will  not  be 

o 

constant  with  cr  , of  course,  if  a generator  of 
finite  impedance  is  taken  into  account. 

Errors  due  to  the  assumption  of  zero  frequency  may 
be  expected  to  become  appreciable  when  (a)  the  series 
inductance  of  the  electrodes  begins  to  play  an  important 
part  in  the  propagation  constant  of  the  transmission  line 
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(the  series  inductive  impedance  becomes  equal  to  the 
series  resistance  at  about  12  c/s),  (b)  the  skin 

depth  of  wave  propagation  in  the  conducting  medium  becomes 
comparable  with  the  dimensions  of  the  system  (the 
frequency  for  a skin  depth  of  10  meters  is  840  c/s, 
that  for  a skin  depth  of  200  meters  is  2.1  c/s). 

Which  of  these  rough  criteria  would  be  the  best  index 
to  the  different  behavior  at  higher  frequencies  must 
await  a more  detailed  investigation  of  the  A.  C.  case. 

Yk  i >L 

M.  S.  Steinberg 
R.  fy.  Jackson 
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